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1. Introduction 

In a recent paper Marques et al[l] derived two new simple supersymmetric partner potentials and showed that the 
ground-state eigenfunction of one of the Hamiltonian operators can be obtained exactly. In order to calculate the 
remaining states they resorted to the Rayleigh-Ritz variational method. However, the convergence rate of their 
approach appears to be insufficient to reveal the relation between the eigenvalues of the partner Hamiltonians 
clearly. 

The authors also estimated the ground-state energy of the quartic anharmonic oscillator by means of perturbation 
theory starting from the known solution for one of the partner Hamiltonian operators. 

The purpose of this paper is twofold. First we improve the variational method proposed by Marques et al[l] and, 
second, we show that the Riccati-Pade method (RPM)[2-4] provides highly accurate results for the eigenvalues 
of those partner Hamiltonian operators. In addition, we show that exactly the same Rayleigh-Ritz variational 
method applied to the SUSY partner Hamiltonians is also suitable for the calculation of the eigenvalues of the 
quartic anharmonic oscillator. 
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In Sec. 2 we apply the Rayleigh-Ritz variational method and the RPM to the pair of SUSY Hamiltonians proposed 
by Marques et al[l] and also to the quartic oscillator. In Sec. 3 we comment on those approaches and draw 
conclusions. 

2. SUSY partners 

The partner Hamiltonian operators are given by[l] 

H± = -£ s +x 4 ±2\x\ (1) 

where we have chosen g = 1 without loss of generality. 
In order to solve the Schrodinger equation 

H ± ^ = Et^t (2) 
approximately, Marques et al[l] chose the nonorthogonal basis set 

tp j (x) = x j - 1 e-M 3/s ,j = l,2,... (3) 

and derived a single secular equation for both the even and odd states. Since states of different parity do not mix 
the dimension of their secular equation is twice the size of what would be actually necessary for a calculation of 
the same accuracy. 

In this paper we propose to work on the half line x £ [0, oo) taking into account the appropriate boundary 
conditions at the origin: ip(0) 7^ 0, i/»'(0) = for the even states and ip(0) = 0, ip'(0) 7^ for the odd ones. The 
appropriate basis sets are 

Pi = x 3 e~ x /3 , j = 0,2,3,... 

P°j = 1.2.3.... (4) 

for the former and latter case, respectively. 

Tables 1 and 2 show Rayleigh-Ritz results for the first eigenvalues of H- and H+, respectively, in terms of the 
size N of the basis sets. We appreciate that the variational results decrease with the number of basis functions as 
expected. However, we are not aware of any proof that either of the basis sets (3) or (4) is complete. We see that 
present results are smaller, and therefore more accurate, than those of Marques et al[l]. In addition, the results 
in their tables 1 and 2 exhibit unchanged entries when the size of the basis set increases from mtom+1 which 
reflects the fact that the added basis function improves the even (odd) state but has no effect on the odd (even) 
one. As a result their matrices are twice the size of what is necessary for a calculation of the same accuracy. 
The results of Tables 1 and 2, as well as those of Marques et al[l], suggest that E~ = E^_ 1 , n = 1,2,..., as 
expected for a SUSY pair of partner Hamiltonian operators. However, those variational results are not sufficiently 
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accurate and, as argued above, we are not aware of any proof of convergence. In order to obtain more accurate 
results we resort to the RPM[2-4]. To this end we consider the modified logarithmic derivative of the wavefunction 

"■>-S-$$ < 6 > 

where s — or s — 1 for even or odd states, respectively. It satisfies the Riccati equation 

f(x) = f(xf - ^f(x) + E~x 4 T 2x,x>0 (6) 
and can be expanded in a Taylor series about the origin 

oo 

/(x) = X>^ (7) 

3=0 

where the coefficients fj can be easily obtained from the Riccati equation (6). We have fo = —i/)'(p)/ij)(0) = for 
the even states and choose fo = for the odd ones in order to remove the unbalanced term 2fo/x that appears 
when we substitute the expansion (7) into the Riccati equation (6). Note that this application of the RPM to even 
potentials is different from the one in earlier papers [2-4] because in the present case we restrict the calculation 
to the half line. The reason for such a modification is that the present even potentials cannot be expanded in an 
a; 2 -power series. 

With the coefficients of the expansion (7) we construct the Hankel determinants Hp = \fj+j+d-\ |fj = i that are 
polynomial functions of the energy and obtain the approximate eigenvalues from the roots of H%(E) = 0[2-4]. 
More precisely, we expect that sequences of roots El^'^ converge towards the eigenvalues E n , n — 0, 1 , . . ., of 
the corresponding Hamiltonian operator as D increases. Since the SUSY partner potentials on the half line 
V±(x) = x 4 ± 2x satisfy V+(—x) = ^-(a:) then the RPM yields the eigenvalues of both partner Hamiltonians 
simultaneously [5]. More precisely, the sequences of roots for the even and odd wavefunctions will converge towards 
Eq <£■(!< E% < E^ < . . . and E± < < < Eg < respectively. 

Table 3 shows the lowest eigenvalues for the SUSY partner Hamiltonians estimated from Hankel sequences with 
D < 30, d = and d = 1. The RPM yields the ground-state eigenvalue of H- exactly (Eq = 0) and there is no 
doubt from present accurate eigenvalues that E~ = E^_ 1 . 

Marques et al[l] estimated the ground-state eigenvalues for the Schrodinger equation with the potentials V(x) = x 4 
and V+ (x) from the exact ground-state solution ip^ , Eg for H- . They resorted to logarithmic perturbation theory 
but their results were not accurate. Here, we point out that they could have obtained more accurate results by 
means of the variational method discussed above. In Table 4 we show the results coming from the more efficient 
Rayleigh-Ritz variational method proposed in this paper and we also add the accurate eigenvalues produced by 
the RPM for comparison. 
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3. Conclusions 

It is well known that the Rayleigh-Ritz variational method is more efficient if one treats every symmetry species 
separately. In the present case we have carried out the separation into even and odd functions on the half line by 
means of appropriate boundary conditions at origin. Thus, the size of present secular equations, and consequently 
the computational effort for a given accuracy, is one half that of Marques et al[l]. 

Working on the half line and taking into account the boundary conditions at origin also enables one to obtain 
highly accurate eigenvalues by means of the RPM. In the present case one has to be cautious because the same 
Hankel determinant exhibits roots close to the eigenvalues of both partner Hamiltonian operators. However, it 
is not difficult to identify the sequence of roots that converges to a chosen eigenvalue. The present RPM results 
confirm that the two Hamiltonian operators proposed by Marques et al[l] are already SUSY partners. 
The RPM was applied to sequences of SUSY partner potentials in the past [6]. However, the main interest in that 
paper was to obtain excited-state eigenvalues with Hankel determinants of relatively small dimension. Besides, 
those potentials can be expanded in x 2 -power series and can therefore be treated in the usual way. 

References 



1] 


F. 


Marques, 0. Negrini, and A. J. da Silva, J. Phys. A 45, 115307 (2012). 


2] 


F. 


M. 


Fernandez, Q. Ma, and R. H. Tipping, Phys. Rev. A 39, 1605 (1989). 


3] 


F. 


M. 


Fernandez, Q. Ma, and R. H. Tipping, Phys. Rev. A 40, 6149 (1989). 


4] 


F. 


M. 


Fernandez, arXiv:0807.0655v2 [math-ph] 


5] 


F. 


M. 


Fernandez, J. Phys. A 29, 3167 (1996). 


6] 


F. 


M. 


Fernandez, Q. Ma, D. J. DeSmet, and R. H. Tipping, Can. J. Phys. 67, 931 (1989) 



Francisco M. Fernandez 



Table 1. Eigenvalues E from the R.ayleigh-R.itz method 



N n = 



1 



1.970246841 
1.970246137 
1.969640134 
1.969515816 
1.969507628 
1.969507538 



5.765408776 
5.511879703 
5.507908922 
5.507440820 
5.507202146 
5.507178381 



9.488542554 
9.488418664 
9.406978612 
9.394523525 
9.394324659 
9.394287127 



14.13007837 
13.89266195 
13.87435060 
13.85957259 
13.85838650 



7 



19.48962926 
19.02107688 
18.66009505 
18.64895855 
18.64704858 



24.09989393 33.11107172 

24.07406802 29.81927599 35.91298387 

23.84777162 29.26028194 35.65656448 

23.80790714 29.25748454 34.98577701 



Table 2. Eigenvalues from the R.ayleigh-R.itz method 



N 



II - 







1 



2 1.970246841 

3 1.970246137 

4 1.969640134 

5 1.969515816 

6 1.969507628 

7 1.969507538 



5.529040685 
5.514418950 
5.510107538 
5.507493533 
5.507185747 
5.507178915 



9.488542554 
9.488418664 
9.406978612 
9.394523525 
9.394324659 
9.394287127 



14.35721210 
14.00103534 
13.92544875 
13.86389218 
13.85851252 
13.85851126 



n - 



n 



n 



It : 



7 



19.48962926 24.52838046 

19.02107688 24.52809061 33.11107172 37.63335229 

18.66009505 23.91353748 29.81927599 36.33167441 

18.64895855 23.81074205 29.26028194 35.02577244 

18.64704858 23.81074194 29.25748454 34.98424967 



Table 3. Eigenvalues for the partner SUSY Hamiltonian operators from the R.iccati-Pade method 



Even = 0) 



Odd (s = 1) 





1.9695075137502948249 
5.5071777771459699676 
9.3942674378738914658 
13.858371936541300147 
18.645975633988444799 
23.8071859179857669 
29.23255455062149608 
34.94633571889061 



1.9695075137502948249 
5.5071777771459699676 
9.3942674378738914658 
13.858371936541300147 
18.645975633988444799 
23.807185917985766918 
29.2325545506214961 
34.9463357188906106 
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Table 4. Eigenvalues E n for the quartic oscillator 



N 


n = 


n = 1 


n = 2 


2 


1.077335422 


3.804924324 


8.102531212 


3 


1.061889825 


3.802939362 


7.486487866 


4 


1.060417725 


3.800572274 


7.456738018 


5 


1.060362727 


3.799746874 


7.456540880 


6 


1.060362727 


3.799674368 


7.455826496 


7 


1.060362223 


3.799673299 


7.455703721 


RPM 


1.0603620904841828996 


3.7996730298013941688 


7.4556979379867383922 




n = 3 


n = 4 


n = 5 


2 


11.86759677 






3 


11.77079816 


16.91718138 


21.92232689 


1 


11.67894913 


16.29307103 


21.79883488 


5 


11.64634787 


16.29029658 


21.28503933 


6 


11.64483615 


16.26577116 


21.24138553 


7 


11.64480669 


16.26191320 


21.24054440 


RPM 


11.644745511378162021 


16.261826018850225938 


21.238372918235940024 



